Source of Acquisition 
NASA Goddard Space Flight Center 


The 2004 Hyperflare from SGR 1806-20: Further Evidence for 

Global Torsional Vibrations 

Tod E. Strohmayer 1 and Anna L. Watts 2 

ABSTRACT 

We report an analysis of the archival Rossi X-ray Timing Explorer (RXTE) 
data from the December 2004 hyperflare from SGR 1806-20. In addition to 
the « 90 Hz QPO first discovered by Israel et al., we report the detection of 
higher frequency oscillations at » 150, 625, and 1,835 Hz. In addition to these 
frequencies there are indications of oscillations at « 720, and 2,384 Hz, but with 
lower significances. The 150 Hz QPO has a width (FWHM) of about 17 Hz, 
an average amplitude (rms) of 6.5%, and is detected in average power spectra 
centered on the rotational phase of the strongest peak in the pulse profile. This is 
approximately half a rotational cycle from the phase at which the 90 Hz QPO is 
strongly detected. The 625 Hz oscillation was first detected in an average power 
spectrum from nine successive cycles beginning approximately 180 s after the 
initial hard spike. It has a width (FWHM) of ~ 2 Hz and an average amplitude 
(rms) during this interval of 9%. We find a strong detection of the 625 Hz 
oscillation in a pair of successive rotation cycles begining about 230 s after the 
start of the flare. In these cycles we also detect the 1,835 Hz QPO with the 
625 Hz oscillation. The rotational phase in which the 625 Hz QPO is detected is 
similar to that for the 90 Hz QPO, indeed, this feature is seen in the same average 
power spectrum. During the time the 625 Hz QPO is detected we also confirm 
the simultaneous presence of 30 and 92 Hz QPOs, first reported by Israel et al. 

The centroid frequency of the 625 Hz QPO detected with RXTE is within 1 Hz 
of the « 626 Hz oscillation recently found in RHESSI data from this hyperflare 
by Watts & Strohmayer, however, the two detections were made in different 
phase and energy intervals. Nevertheless, we argue that the two results likely 
represent detections of the same oscillation frequency intrinsic to the source, but 
we comment on some of the difficulties in making direct comparisons between the 
RXTE and RHESSI measurements. We argue that these new findings provide 
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further evidence for a connection of these oscillations with global seismic modes 
of neutron stars, in particular, the high frequency signals appear to be consistent 
with toroidal modes with n > 0. We discuss their implications in the context of 
this model, in particular for the depth of neutron star crusts. 

Subject headings: stars: magnetic — -pulsar: individual (SGR 1806-20) — stars: 
neutron— stars: rotation— stars: oscillations — X-rays: stars 

1. Introduction 

Recent studies of high time resolution data from two magnetar hyperflares recorded by 
the Proportional Counter Array (PCA) onboard the Rossi X-ray Timing Explorer (RXTE) 
have resulted in the discovery of a new phenomenon associated with these events. Both 
the December, 2004 event from SGR 1806-20, and the August, 1998 hyperflare from SGR 
1900+14 produced fast, rotation-phase-dependent X-ray oscillations in the 20 - 150 Hz range. 
First, Israel et al. (2005) reported the discovery of ~ 18, 30, and 90 Hz quasi-periodic 
oscillations (QPO) in the December, 2004 event, and suggested that the 30 and 90 Hz QPOs 
could be linked with seismic vibrations of the neutron star crust. Strohmayer & Watts (2005) 
then reported the discovery of a sequence of QPOs in the SGR 1900+14 hyperflare. They 
found a set of frequencies; 28, 53.5, 84, and 155 Hz, which could plausibly be associated 
with a sequence of low l toroidal modes (denoted rfo) of the elastic neutron star crust (see, 
for example, Hansen k, Cioffi 1980; McDermott, van Horn & Hansen 1988; Duncan 1998; 
Piro 2005). In both hyperflares the oscillations are episodic, that is, their amplitudes vary 
considerably with time and rotational phase. Most recently, Watts &; Strohmayer (2006) 
examined Ramaty High Energy Solar Spectroscopic Imager (RHESSI) data for the SGR 
1806-20 event, confirmed the presence of the 18 Hz and 90 Hz QPOs, and found evidence 
for additional oscillations at 26 Hz and 626 Hz. 

The similar phenomenology of the oscillations in the two sources, as well as the closeness 
of some of the measured frequencies argues rather convincingly that we are seeing the same 
physical process in each case. The connection with torsional modes of the crust seems 
plausible for several reasons; 1) The observed frequencies are consistent with theoretical 
expectations for such modes, and can be more or less self-consistently associated with a 
sequence of modes with varying spherical harmonic index, l. 2) The magnetic instability 
which powers the hyperflares is likely associated with large scale fracturing of the neutron 
star crust (Flowers & Ruderman 1977; Thompson & Duncan 1995; Duncan 1998; Thompson 
& Duncan 2001; Schwartz et al. 2005), and will almost certainly generate seismic motions 
within the star. 3) The strong phase dependence argues for a mechanism associated with 
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particular sites on the stellar surface, such as a fracture zone or magnetic field bundle. 
4) Mechanical motions provide a natural explanation for the relatively high coherence of 
the oscillations. 5) Horizontal motions of the crust could modulate the beaming pattern 
associated with the strong magnetic field, providing a mechanism to modulate the X-ray flux. 
Moreover, beaming can act as an “amplifier,” producing potentially large X-ray modulations 
from modest horizontal displacements. Although the present evidence for torsional modes 
is very suggestive, it is not yet definitive. Levin (2006), for example, argues that toroidal 
modes may damp too quickly to account for the detection of oscillations some minutes after 
the onset of the flare. In addition some of the detected frequencies (18 Hz, 26 Hz) do not 
fit easily into current torsional mode models, without invoking magnetic splitting or other 
complications. 

In our previous study of the RXTE data from the SGR 1900+14 hyperflare we computed 
average power spectra for different rotational phases and time intervals during the event. 
This method should be more sensitive to oscillations which are localized in rotational phase, 
as borne out in the SGR 1900+14 study. In this paper we report the results of a similar, 
phase averaged timing study on the now public archival RXTE data from the SGR 1806- 
20 hyperflare. Using this method we find strong evidence for additional, higher frequency 
oscillations during the flare. 


2. Observations and Data Analysis 

An overview of the RXTE data recorded from the December, 2004 hyperflare from SGR 
1806-20 is given by Israel et al. (2005). Data were recorded in the “Goodxenon_2s” mode 
that allows for time resolution up to « 1 fxs. The hyperflare intensity profile, folded at the 
rotational frequency (see Figure 1), shows 3 peaks. In subsequent discussions we will refer 
to these as peaks 1, 2, and 3, in order of decreasing peak intensity, respectively. We will also 
refer to the region of pulse phase between peaks 2 and 3 as the “interpulse” region. 


2.1. Kilohertz Oscillations 

Since previous studies have shown the « 90 Hz QPO to be extremely robust, and its 
approximate location in rotational phase has been linked to the interpulse region, by both 
Israel et al. (2005) and Watts & Strohmayer (2006), we began our study by computing 
average power spectra around this phase range, but for a sequence of different time intervals 
during the flare. We used 3 s intervals to compute each individual power spectrum and we 
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initially restricted our search to a Nyquist frequency of 4096 Hz, although these PC A data 
can, in principle, be used to sample to much higher frequencies. We found a significant signal 
at 625 Hz in power spectra computed from approximately the last third of the hyperflare. 
We now discuss the detection, significance and properties of this oscillation. 

Figure 2 shows an average power spectrum computed from 9 successive rotations starting 
approximately 190 s after the onset of the hyperflare. We used 3 s intervals to compute the 
power spectra. Both the time interval during the flare (top), and the phase region (bottom, 
dashed lines) used in computing this spectrum are shown in Figure 2. This power spectrum 
has a prominent feature at 625 Hz, and the QPO at 92 Hz is also clearly detected. To 
estimate the significance of the 625 Hz feature we extended the frequency range to 65 kHz, 
and used all frequency bins above 800 Hz to estimate the noise power distribution. Figure 
1 shows the power spectrum to 65 kHz (main panel), as well as the distribution of noise 
powers (inset panel). We fit the noise power histogram with a % 2 function, and found a 
reasonable fit for 87 degrees of freedom (dof), and a small reduction in the Poisson level of 
0.01. This function is also plotted in the inset panel of Figure 2 (solid curve). Using this 
noise distribution we find a single trial significance of 7 x 10 -11 for the 625 Hz feature. The 
power spectrum has 1,536 bins up to 4096 Hz (2.66 Hz resolution), which was the top end 
of our search range. We emphasize that we extended the frequency range only after the 
search, simply to better characterize the noise power distribution. This gives a probability 
of 1.1 x 10 -7 to find a peak this high in the power spectrum. The spectrum shown in Figure 
1 was one of a sequence from across the entire duration of the flare. We computed these by 
overlapping the time intervals, so each spectrum in the sequence is not fully independent. 
However, even using the total number of spectra computed (47), we have a significance 
< 1 x 10~ 5 , so this is a robust detection. We fit the QPO with a Lorentzian profile and find 
a center frequency and width of 625.5 ± 0.15 Hz, and 1.8 ± 0.4 Hz, respectively. The 625 
Hz QPO has an average amplitude during this interval of 8.5%, and a width (FWHM) of 
about 2 Hz. We further note that the strength of the 625 Hz feature is comparable to that 
of the 92 Hz signal, which provides additional confidence that it is not simply a statistical 
fluctuation. 

We searched the RHESSI data to see if there was any indication of a simultaneous 625 
Hz peak but did not detect anything significant at this time. This is however not surprising 
given the amplitude, and the fact that for such high frequencies we can only use the RHESSI 
front segments, which have a much lower countrate than RXTE (see the discussion in Watts 
& Strohmayer (2006)). In addition the RXTE signal is strongest at energies less than 45 
keV; at low energies the RHESSI spectrum exhibits much higher background levels, which 
could swamp such a weak signal. 
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We computed dynamic power spectra in the vicinity of 625 Hz to search for time de- 
pendence of the signal. The results suggest that the 625 Hz oscillation is most strongly 
associated with the falling edge of peak, 2 of the pulse profile, and that it’s amplitude can 
vary strongly with time. For example, Figure 3 shows a portion of the hyperflare that con- 
tributes strongly to the 625 Hz signal detected in the average power spectrum described 
above. Contours of constant Fourier power in the vicinity of 625 Hz are plotted over the hy- 
perflare profile. We used 1/2 s intervals and started a new interval every 0.1 s. One pulsation 
cycle is plotted. The strong correlation of power contours with the falling edge of peak 2 is 
rather similar to the behavior seen in the 84 Hz oscillation from SGR 1900+14 (Strohmayer 
& Watts 2005, see their Figure 3), and these similarities suggest it may be connected with 
the X-ray modulation mechanism. 

We also computed several average power spectra in the immediate vicinity of the strong 
625 Hz signal seen in the dynamic power spectrum. An average from two consecutive cycles 
beginning with the cycle prior to that shown in Figure 3 reveals evidence for an oscillation 
at an even higher frequency of 1,840 Hz. The average power spectrum from these two cycles 
is shown in Figure 4. The two spectra are the same except for being plotted at different 
resolutions. The peak at 1,840 Hz is most prominent in the top trace, at 5.33 Hz resolution 
(the bottom trace has 2.66 Hz resolution). Both the 92 Hz and 625 Hz peaks also stand 
out prominently in this spectrum. We estimated the significance in the manner discussed 
above, and found chance probabilities of 1.1 x 10~ 6 , and 1.1 x 10 -5 for the 625 and 1,840 Hz 
peaks, respectively. This is an estimate of the probability of finding a peak as high as those 
observed in this power spectrum up to 4096 Hz. We computed eight such spectra, so this 
will reduce the significances somewhat, but again, the detections are further supported by 
the comparable amplitude of these peaks to that of the well established 92 Hz oscillation, 
which is also detected in the same power spectrum. 

We also see weaker indications for a few additional high frequency peaks during this 
time interval. Figure 5 compares two power spectra, one (top) from the two cycle interval in 
which the 1,840 Hz signal is seen, and the other (bottom) an average of two cycles beginning 
six cycles earlier than the cycle shown in Figure 3. The frequency resolution is 10.66 Hz in 
both plots. The top spectrum is actually the same as in Figure 4, but plotted with lower 
resolution. The 625 and 1,840 Hz peaks still stand out prominently. We estimated the 
standard deviation in the noise distribution as before, and we placed two horizontal dotted 
lines 4a (single trial) above the mean in each case. In the bottom spectrum, there are also 
two peaks which reach above the 4cr line. These have frequencies of 720 and 2,384 Hz. 
We estimate significances (single trial) of 4.4 x 10 -6 and 1.3 x 10 — 5 , respectively, for these 
peaks. A third peak, at 976 Hz, does not reach the 4er line, but is a 3cr deviation, with an 
estimated significance of 7.6 x 10 -4 . While we do not claim that all these peaks represent 
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strong detections, it does seem unlikely by chance to have several such sizable deviations 
from the mean in a single spectrum, and thus we think this could be an indication of the 
presence of additional, weaker signals in the data. 

We also searched the RXTE data for evidence of a 626 Hz QPO earlier in the tail, when 
a QPO of that frequency was seen in the RHESSI data (Watts & Strohmayer 2006). We 
were not able to find a corresponding signal in the PC A data. However, this is not entirely 
surprising. The ■ 626 Hz signal was only observed in the RHESSI data for photons with 
recorded energies in the range 100-200 keV. Whilst RHESSI is sensitive at these energies, 
the PCA does not record photons with energies above « 100 keV. Comparing spectra from 
the two spacecraft is difficult because of the scattering caused by the oblique incidence angles. 
However we can get some idea of the magnitude of the energy shift between the two datasets 
by comparing the positions of the spectral maxima. For the PCA the maximum is at 20 keV, 
for RHESSI it is at 40 keV. If incident photons are downshifted in energy by only 20keV for 
the RXTE PCA compared to RHESSI, then we would not expect to see the RHESSI 626 
Hz signal, since the PCA does not capture sufficient high energy photons. Unfortunately 
the RXTE High Energy X-ray Timing Experiment (HEXTE), which covers a higher energy 
range than the PCA, recorded the main flare but not the decaying tail (Smith et al. 2005). 


2,2, 150 Hz Oscillation 

We also computed average power spectra centered on the rotational phase of peak 1. 
Figure 2 shows the range of phases used (dashed-dot lines). Average power spectra centered 
on this phase over the entire hyperflare show a significant QPO centered near 150 Hz. Figure 
6 shows the power spectrum in which we detected this QPO. We averaged 45, 3 second power 
spectra, and the frequency resolution is 8 Hz. The peak near 150 Hz is evident. We modelled 
the power continuum in the range above 30 Hz with a constant plus a power law, and then 
rescaled the spectrum by dividing by the best fitting continuum model. We then estimated 
the noise power distribution in the same way as described earlier, by fitting a % 2 distribution 
to a histogram of noise powers in the rescaled spectrum. Based on this we estimate a 
significance for the highest bin in the QPO profile as 2.8 x 10~ 5 , which includes the number 
of frequency bins searched (up to 4096 Hz). This is a very conservative estimate of the 
QPO significance, because it only uses the highest bin, and the QPO peak is resolved at this 
frequency resolution. We fit the QPO profile with a lorentzian and find a centroid frequency 
of 150.3 ± 1.6 Hz, a width of 17±5 Hz, and an amplitude (rms) of 6.8 ±1.3%. If we construct 
averages from the first and second halves of the hyperflare, we still detect the 150 QPO, so 
this oscillation seems to persist through most of the hyperflare. This feature is not detected 
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in the RHESSI data, but this is not unexpected given the low amplitude of the QPO and 
the lower countrates of the RHESSI front segments. 

2.3. Lower frequency QPOs 

Israel et al. (2005) reported the marginal detection of 18 Hz and 30 Hz QPOs over a 
period late on in the burst , without exploring for any phase dependence. Watts & Strohmayer 
(2006) then found evidence in the RHESSI dataset for broad QPOs at 18 and 26 Hz, with 
a weaker feature at 30 Hz. These QPOs showed a strong rotational phase dependence, 
appearing at the same phase as the 92.5 Hz QPO. For these lower frequency QPOs we 
were able to include photons detected by the RHESSI rear segments, so that the RHESSI 
countrate exceeded that recorded by RXTE. 

The broad QPOs at 18 Hz and 26 Hz were seen in the RHESSI dataset between 50 and 
200s after the main flare. We were interested to see whether these QPOs were detectable 
in the RXTE data. Folding up exactly the same time period and rotational phase used 
in the RHESSI analysis we were unable to make a statistically significant detection, not 
unexpected given the discrepancy in countrates and the low fractional amplitude of the 
QPOs. We therefore focused on two shorter periods in the RHESSI data when the 18 Hz 
and 26 Hz oscillations were particularly strong, and searched for them at corresponding times 
in the RXTE data. In both of these cases corresponding strong peaks are observed in the 
RXTE data, confirming the detection of QPOs at these frequencies. Interestingly the 26 Hz 
QPO was stronger in this period in the RXTE dataset than it is in RHESSI, suggesting that 
other factors may offset the higher countrate (see below). 

What about the 30 Hz QPO? If we consider the same time interval indicated in Figure 
1, but shift the phase range over which we compute power spectra to later phases, that is, 
exclude the falling edge of peak 2, the 625 Hz signal drops to undetectable levels. Figure 7 
shows the average power spectrum of this “interpulse” region during this time interval, and 
confirms the initial claims of Israel et al. (2005) that signals at « 30 and 92 Hz are present 
towards the end of the hyperflare. These QPOs are extremely significant, for example, we 
conservatively estimate the significance of the 30 Hz QPO at 1.4 x lO -11 , and the 90 Hz 
feature is even more significant. These QPOs can be fitted with lorentzian profiles. Based 
on such fits we find centroid frequencies of 28.98 ± 0.4, and 92.9 ± 0.2 Hz; quality factors 
(uq/Au) of 7±0.8, and 39 ±5; and average amplitudes of 20.5 ±3, and 19.2 ±2, respectively. 

The 30 Hz feature in the RXTE spectrum is far more significant than the weaker feature 
at 30 Hz found in the RHESSI spectrum. This is at first glance surprising, since RHESSI has 
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a higher countrate once the rear segments are included. There are however several factors 
that probably contribute to this discrepancy. Firstly, we are in the interpulse region, where 
RHESSI’s high background levels are more important. Secondly, one cycle of this time period 
was excised from the RHESSI data because of an artificial spike in countrate as a protective 
attenuator was removed. In addition, as discussed above, we find occasions where the 26 Hz 
QPO appears far more prominently in the RXTE data than in the RHESSI data, suggesting 
that other factors such as scattering off the spacecraft and the Earth may reduce RHESSI’s 
effectiveness and offset the higher countrate. 


3. Discussion 


Using phase averaging techniques we have found evidence for additional, higher fre- 
quency QPOs during the December, 2004 hyperflare from SGR 1806-20. These additional 
frequencies and the time and phase dependence of the detections suggest an extremely rich 
and complex phenomenology is associated with the magnetar oscillations. 

The detection of a strong 625 Hz QPO in the RXTE data is particularly interesting. 
The frequency is remarkably similar to that of a QPO found by us in the RHESSI data, 
but there are some intriguing differences in the properties of the two QPOs. In addition to 
having a slightly lower centroid frequency (there is frequency overlap of the QPO profiles 
however), the RXTE QPO has substantially lower fractional amplitude and coherence value 
Q. The RXTE QPO appears later in the tail of the flare, and the apparent photon energies 
of the two QPOs also differ (see Figure 8). 

The RHESSI detection used events from the front segments, in the 100 - 2100 keV energy 
band (see Watts & Strohmayer (2006) for a discussion of RHESSI data in this context). 
The phase interval was centered on the rising edge of peak 1 (Watts & Strohmayer (2006), 
Figure 4), and the average was obtained from about 50 to 200 s after the start of the 
hyperflare. An average power spectrum from RXTE for the same time and phase ranges 
shows no corresponding signal in the PCA. However, it is difficult to know exactly how to 
compare the energies of events detected in the two instruments. For RHESSI, many of the 
events observed in the front segments were direct (ie. meaning that they were not scattered 
before detection), wheareas this is likely not the case with RXTE, as virtually no photons 
would have made it directly into the PCA xenon volumes without having scattered off of 
some part of the spacecraft. It is quite possible that no photons with incident energies in the 
100 - 200 keV range were detected by the PCA, which would account for the non-detection. 
The later, lower amplitude RXTE QPO is not detected in the RHESSI data, most probably 
due to the lower countrates. 
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The lack of simultaneous detections is therefore disappointing, but not surprising. Per- 
haps more important is the fact that the two instruments detected signals independently at 
consistent frequencies. The QPO profile of the 625 Hz feature detected with RXTE over- 
laps the 1 Hz bin (centered on 626.5 Hz) in which the signal was detected with RHESSI. 
Although it is difficult to determine an a priori probability, it seems extremely unlikely that 
both instruments would detect signals at a consistent frequency unless that frequency were 
intrinsic to the source. We therefore think there is compelling evidence to associate this 
frequency with SGR 1806-20. 

The finding of additional characteristic frequencies in the SGR 1806-20 hyperflare further 
strengthens the connections between the oscillations seen in both the SGR 1806-20 and SGR 
1900-1-14 events. In each case, oscillations near ~ 30, ~ 90 and « 150 Hz have been seen. 
This strongly suppports the notion that the same physical processes are involved in these 
objects, and makes sense in the context of global oscillation modes; if the neutron stars have 
similar masses and magnetic fields, then they should ring with the same set of characteristic 
frequencies. 

The 30 Hz, 92 Hz and 150 Hz QPO can be interpreted within the global oscillation model 
as n = 0, 1=2, 7 and 12 torsional modes of the neutron star crust. Watts & Strohmayer 
(2006) suggested that the 626 Hz oscillation seen in RHESSI data could be associated with 
a toroidal mode having a single node in the radial displacement eigenfunction (fa modes). 
Recent estimates from Piro (2005, see his Figure 3) suggest that this frequency is consistent 
with theoretical expectations. Our confirmation of the same frequency in the RXTE data, 
provides a firm indication that we are seeing not only modes with a sequence of different 
£- values, but also modes with different n- values. 

It is possible that the 625 Hz QPO seen in the RXTE dataset represents an evolution of 
the QPO seen earlier in the flare in the RHESSI data. However, the mode model also offers a 
possible explanation for the detection of two ~ 625 Hz QPOs with differences in properties. 
Calculations (see Hansen &Cioffi 1980; Piro 2005), indicate that the n = 1 mode frequencies 
are not very sensitive to l. As such we could be seeing two n = 1 modes of different l, excited 
at different times by either aftershocks or mode coupling. The implications of detecting an 
n = 1 mode are discussed in more detail below. First, however, let us consider how the other 
frequencies fit into the current oscillation models. 

The 18 Hz and 26 Hz QPOs are at first glance too low in frequency to fit simple 
torsional mode models, however, mode calculations with realistic magnetic fields will have 
ra-dependent splitting, perhaps producing an observable multiplet structure. This might 
accomodate the observed spread from about 26 - 30 Hz in a putative 2 to mode, but detailed 
Calculations will be required to explore this quantitatively. Global magneto-elastic modes 
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may also be a possibility here (see for example Glampedakis, Samuelsson & Andersson 
(2006)), although there is an issue of damping (see below). 

The stronger high frequency QPOs at 625 and 1840 Hz can be nicely accomodated 
within the torsional mode hypothesis as low l, n = 1 and n — 3 modes, respectively (see, for 
example, Figure 3 in Piro 2005). The less significant features at 720, 976, and 2,384 Hz, may 
also fit within this scenario. The 720 Hz might represent a low n mode with significantly 
different l , as pointed out by Hansen & Cioffi (1980), the l dependence of such modes is not 
strong, but can produce a scatter in frequency on the order of 15%. The other frequencies 
might correspond to modes with n = 2, and n > 3, respectively. In addition, at these 
high frequencies, other modes, such as the crustal spheroidal modes (see McDermott et al. 
1988), may become relevant. Again, more detailed mode calculations, including the effects 
of magnetic fields, will be needed to make more quantitative mode identifications. 


3.1. Mode Excitation and Damping 

In a recent paper Levin (2006) has argued that the high frequency magnetar oscillations 
cannot be due to mechanical modes of the crust. The crux of Levin’s argument is that such 
modes would damp very quickly (< 1 s) via Alfven radiation into the neutron star core. 
While the issue of damping is of course an important one we think the conclusion that crust 
modes cannot be driving the oscillations is very premature, especially in fight of the large gaps 
in our understanding of the magnetic field structure inside neutron stars. For example, the 
damping into the core will depend sensitively on the amplitude at the crust core boundary, 
and previous calculations show that this amplitude can be substantially smaller than that 
at the top of the crust (see McDermott, van Horn & Hansen 1988, for example). Thus, 
approximating the amplitude as a single constant for a given mode may not be very accurate. 
The physics at the boundary layer at the base of the crust will be crucial (see for example 
Kinney & Mendell (2003)). Moreover, the damping will likely depend rather sensitively 
on both the magnetic field geometry and the particular displacement pattern of individual 
modes, so some modes might be much more strongly damped than others. Additionally, 
the current observational evidence suggests that the mechanism which modulates the X- 
ray flux is strongly “tunable.” That is, at times when certain frequencies are detected their 
instantaneous amplitudes can be quite high, almost certainly higher than the relative surface 
motions. So, it may be that the modulation mechanism can amplify relatively small surface 
motions. Also, if mode energy is radiated into the core it may be possible for some of this 
energy to be pumped back into mechanical motions of the crust, producing a kind of feedback 
effect (Duncan 1998). Thus, it may be that the excitation of modes occurs over a time scale 



not too dissimilar from that of the hyperflare itself. 


3.2. Implications for Neutron Star Crusts 

If the 625 Hz QPOs are indeed n = 1 torsional modes, the implications for the study 
of neutron star structure are profound. The ratio of frequencies of the n = 1 and n — 0 
torsional modes enables one to deduce the depth of the neutron star crust (see Hansen & 
Cioffi 1980; Piro 2005). The basic physics behind this estimate is that the n — 0 mode 
frequencies are set by the horizontal wavelength, which scales with the stellar radius, R, 
whereas the n > 0 mode frequencies are dominated by the vertical wavelength, which senses 
the depth of the crust, A R. In the limit of a thin crust and a constant shear wave speed, one 
finds the simple expression (see Hansen & Cioffi 1980; McDermott et al. 1988; Piro 2005), 

(i(i+l)) 1/a AJi 
^ti>o 3n R 

If we identify 30 Hz with the 1 = 2, n = 0 mode frequency, and 625 Hz with the n = 1 
mode (and l « n), the above expression predicts A R/R = 0.06. In most fnodern neutron 
star models, however, the crust is not “thin” in the sense of A R/R « 1, and the shear 
speed has some variation with depth, so this simple estimate should probably be considered 
a lower limit. A more accurate estimate can be obtained by comparisons with more detailed 
mode calculations. For example, using the calculations of McDermott et al. (1988) we 
evaluated the expression, (i/i^AR) / (i/ 2 t oR) , that according to the approximate theory should 
be the proportionality constant between the mode frequency ratio and the fractional depth 
of the crust. For the four models tabulated by McDermott et al. (1988) we find that this 
expression ranges from 2.15 to 2.65. Adopting this range as representative, and using the 
scale factors to convert our observed mode frequency ratios into depths, we obtain the range, 
0.1 < A R/R < 0.127. We note that the models from McDermott et al. use softer equation 
of state models than are currently favored, so it is important to repeat this exercise with 
more up to date neutron star models, and also incorporate magnetic fields. 

Measuring the thickness of a neutron star crust also conveys information on the equation 
of state (EOS). In other words, stellar models of a given mass computed with different EOSs 
will in general have crusts of different depths. For example, Lattimer (2006) has recently 
outlined how constraints on the fractional depth and stellar mass could be used to constrain 
the nuclear symmetry energy and the nuclear force model. The crust depth dependence 
on the EOS also has implications for pulsar glitch models (see, for exmaple, Crawford & 
Demianski 2003; Link, Espstein & Lattimer 1999). 
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Neutron star seismology, as revealed by magnetar hyperflares, has great potential as a 
probe of stellar structure, composition and magnetic field geometry. It is clear however that 
serious theoretical development is required - the coupling of the modes to the magnetic field, 
and the modulation of the lightcurve being two of the most pressing areas of investigation. 
The complex temporal variation, however, also requires explanation. 

Figure 8 shows the periods during the tail of the flare when the different QPOs are de- 
tected in either the RXTE or RHESSI datasets. Some of the oscillations are seen throughout 
the flare, others are detected only after the rise in emission half way through the tail, and 
the highest frequency modes are very short-lived. The data suggest that some QPOs are 
excited by the main flare, whereas others are excited later in the tail, perhaps by aftershock 
fracturing and reconnection. The association of several of the QPOs with the rise in emission 
late in the tail lends support to the idea of continued excitation. However, there remains the 
possibility that all of the QPOs are excited by the main flare, with their visibility depend- 
ing on beaming effects that vary during the tail. Beaming will be influenced by both the 
magnetic field and by the plasma surrounding the magnetar. Selection by beaming might 
also account for the fact that -only some of the myriad oscillations predicted by theory are 
observed - although mode excitation will also depend sensitively on the speed and extent 
of the fracture. Some modes may also damp far too quickly to be detected, depending on 
local conditions in the crust and magnetic field geometry. We believe that these factors are 
more than capable of accounting for the fact that only a few frequencies reach detectable 
amplitudes, an issue of concern to Levin (2006). 

In the periods when they are detected, the QPOs are far from static. Their amplitudes 
wax and wane, and there is evidence for frequency shifts and possible multiplet splitting. In 
Figure 9 we show a sequence of dynamical power spectra for the Peak 2/Interpulse region, 
for the time period when the 92.5 Hz QPO is active. The presence of the strong QPO at this 
frequency from 180-230s after the main flare is clear. If we look at the cycles around this 
period however, there are hints of interesting behavior. Back at only 80s after the flare, weak 
detections are made, with a frequency that starts at around 78 Hz and rises up to a strong 
92.5 Hz QPO that persists for a couple of cycles from 120-130 s. The QPOs then weaken, 
and the frequency seems to fall again by a few Hz, until the strong 92.5 Hz signal reappears 
at 180s. At later times, as already noted by Israel et al. (2005), the frequency rises to 95 Hz 
and possibly higher, as the amplitude falls and the signal becomes undetectable. There is 
also some evidence for frequency movement in the lower amplitude QPOs, although this is 
harder to detect because the QPOs are weaker and the amplitudes more variable. As such 
it is not possible to verify whether the low QPO frequencies move in concert with the higher 
frequency QPO. 
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The variation of amplitude over the rotational cycle, and the consequent interruption of 
the signals, means that we cannot rigorously connect cycles, as one would wish to in order 
to verify that the frequency is evolving. But our results are suggestive of this eventuality 
and it merits further study. Evolution of the magnetic field in the aftermath of the flare, 
relaxation of the deep crustal structure, and the evolution of the surrounding plasma, are 
all candidates for causing such frequency shifts. Alternatively we could be seeing different 
members of a magnetically-split multiplet excited, perhaps, by mode coupling. 


4. Conclusions 


Our study of the RXTE data from the SGR 1806-20 hyperflare indicates that a complex 
pulsation phenomenology is associated with magnetar hyperflares. The discovery of new kHz- 
range frequencies consistent with theoretical predictions for n > 0 torsional modes provides 
strong evidence that we may be seeing vibration modes of the neutron star crust excited 
by these catastrophic events. If this is true, then it opens up the exciting prospect of 
probing the interiors of neutron stars in a manner analogous to helioseismology. Additional 
excitement is warranted when we consider that all the current datasets used to explore these 
oscillations have been purely serendipitous. That is, they have not been optimized in any 
way for studying these signals. This suggests that a wealth of additional information would 
likely be found from instruments better optimized to capture with high time resolution the 
flood of X-rays produced by these events. Based on the findings presented here we are also 
strongly convinced that more theoretical work is definitely needed to make more accurate 
mode identifications, to better understand the excitation and damping mechanisms of crust 
modes and how they can couple to the X-ray emission, and to make more precise inferences 
on neutron star structure. 
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Pulse Phase 


Fig. 1. — X-ray intensity of the December, 2004 hyperflare as measured by the PGA (top), 
and the average pulse profile (bottom). The curves include all good events detected in 
PCA channels 10 - 200 (nominal energy band from 4 - 90 keV). The main flare takes place 
approximately 4s prior to the zero on the time scale. The time interval in which the 625 Hz 
oscillation was detected is marked by the vertical dashed lines (top). The vertical lines in 
the bottom plot denote the phase ranges used to compute power spectra shown in Figures 
1 (dashed), 5 (dotted), and 7 (dash-dot). The text refers to Peaks 1, 2 and 3: Peak 1 is at 
phase ~ 0.8, Peak 2 at ~ 0.1 and Peak 3 at ~ 0.6. 
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Fig. 2.— Average power spectrum from a portion of the hyperflare from SGR 1806-20 (main 
panel). We averaged nine 3 s power spectra from the time interval marked by the vertical 
dashed lines in Figure 2 (top). The frequency resolution is 2.66 Hz. The inset panel shows 
the distribution of noise powers computed from the frequency range 800 - 65,536 Hz as well 
as the best fitting x 2 distribution (solid). The distribution is plotted in units of o' s. The 
vertical dashed line marks the peak power of the 625 Hz feature. See the text for further 
discussion.' 
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Fig. 3. — Dynamic power spectrum of a portion of the data interval in which the 625 Hz 
oscillation was detected. Contours of constant Fourier power are plotted as a function of 
frequency and time, along with the X-ray intensity as a function of time,. We used 0.5 
s intervals to compute the power spectra, and overlapped the intervals, beginning a new 
one every 0.1 s. A strong contribution to the 625 Hz signal detected in the average power 
spectrum of Figure 1 is evident on the falling edge of peak 2. Contours at Leahy normalized 
power levels of 16, 18, 20, 24, 28 and 32 are shown. A single rotational cycle is plotted. 
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Fig. 4. — Average power spectrum from two consecutive cycles beginning with the cycle 
immediately prior to the pulse displayed in Figure 3 which shows a strong 625 Hz signal. 
Two representations of the same power spectrum are shown, only the frequency resolutions 
differ. The frequency bins are 2.667 Hz wide in the bottom trace and twice that in the top. 
Three frequencies are prominent in the top trace, 92 Hz, 625 Hz, and 1,840 Hz. See the text 
for additional discussion. 
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Fig. 5. — Average power spectra from the time interval in which the 625 Hz oscillation was 
first detected. Shown are average power spectra computed from two consecutive cycles. The 
top panel includes the pulsation cycle shown in Figure 3, and the preceeding cycle (same 
as in Figure 4, only plotted with lower frequency resolution). The bottom panel shows the 
spectrum from earlier in the time interval, beginning six cycles before that shown in Figure 
3. In each case the frequency resolution is 10.66 Hz. In each spectrum a horizontal dotted 
line marks a 4cr (one trial) deviation. The 625 and 1,840 Hz features are detected in the top 
spectrum, while several different frequencies are suggested in the bottom spectrum. See the 
text for more details. 
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Fig. 7. — Average power spectrum from the “interpulse” region during the same time interval 
in which the 625 Hz oscillation was detected. The phase interval used for the average 
spectrum is marked with vertical dotted lines in Figure 2. Strong QPQs at 29 and 92 Hz 
are clearly detected. See the text for additional discussion. 
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Fig. 8. — The time periods when the different QPOs are detectable in either the RXTE 
or RHESSI datasets. Solid lines indicate QPOs that are detected in the Peak 2/Interpulse 
region. Dashed lines indicate QPOs that are detected primarily during Peak 1. The difference 
in properties between the two « 625 Hz QPOs, the earlier one detected by RHESSI and the 
later one by RXTE, is clear. 
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Fig. 9. — The evolution of frequency in the 75-105 Hz band, for the Peak 2/interpulse region. 
Each numbered panel shows a 4s interval for consecutive cycles, the first panel starting « 60 
s after the main flare. The y-axis shows frequency in Hz. Dynamical power spectra, shown 
as contours, are computed from Is segments, overlapping by 0.1s. The minimum contour 
level shown is a power of 15, with increments of 5. Lightcurves are also shown, plotted to 
the same scale on each panel. The plots show amplitude variation, frequency evolution, and 
possible multiplet activity. 




